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Abstract 

Recent CDF data of the inclusive jet cross section shows threshold-like struc- 
tured deviation, around transverse momentum E T (j) ~ 200 ~ 350 GeV. If 
this data is real, not just some statistical fluctuation, is it possible to interpret 
the anomaly in terms of virtual SUSY effects? The purpose of this note is to 
address this question. However, we find that virtual SUSY loop interference 
effects [near the threshold] are too small to explain the CDF data. Our main 
conclusion seems to be on the right track if we assume that the recent global 
analysis of improved parton distributions by Lai et al. is correct. 
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I. INTRODUCTION 



The CDF |I[ and DO collaborations at Tevatron Collider have recently reported 
data for the inclusive jet cross section in pp collisions at y/s = 1800 GeV. Let us recapitulate 
some particulars of this data [we concentrate mainly on the CDF data]. 

• The CDF measurement is based on a data sample of 19.5 pb _1 collected in 1992-93 
with the CDF detector at the Tevatron collider. Jets were reconstructed using a cone 
algorithm. 

• Cosmic rays and accelerator loss backgrounds were removed with cuts on event energy 
timing and on missing transverse energy. The remaining backgrounds are claimed to 
be less 0.5% in any E T bin. 

• The measurements have been reported over a wide range of transverse energy, 15 GeV 
< E T < 440 GeV, and around the central pseudorapidity region 0.1 < \rj\ < 0.7. The 
smallest distance probed is on the order of 10~ 19 m. 

• After accounting for uncertainties the corrected experimental cross section, when com- 
pared to the Next-to-leading [NLO] QCD predictions using MRSD0' parton distribu- 
tion function [PDF's], is significantly higher than the NLO prediction for E T > 200 
GeV. For E T below 200 GeV the agreement between the CDF and the NLO QCD cross 
section is excellent, while the DO results are higher than the NLO PQCD predictions 
within the statistical uncertainties. 

• CDF collaboration have compared their data with other PDF's and a model about 
presence of quark substructure. This will be discussed later. 
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There are basically two logical possibilities for explanations of the CDF data on inclu- 
sive jet production cross sections: (i) the parton distribution functions determined at low Q 2 
region may not be accurate enough to be applicable to the high E T region with E T > 200 
GeV, or (ii) there are some new physics around the electroweak scale. We give a brief review 
of these two possibilities in the following. 

Let us first consider the first possibility. According to the CDF collaboration [||] "the 
excess of data over theory at high E T remains for CTEQ2M, CTEQ2ML, GRV94, MRSA', 
and MRSG parton distribution". The variations in QCD predictions represented a survey of 
then-available distributions. They do not represent the uncertainties associated with data 
used in deriving the PDF's. Inclusion of these new data in a global fit with those from other 
experiments may yield a consistent set of PDF's which will accommodate the high-E^ excess 
within the scope of QCD. Glover et al, [[J conclude that "it is unlikely that the difference 
between the CDF inclusive jet cross section data and standard NLO QCD prediction can be 
attributed to a deficiency in our knowledge of parton distributions". More precisely Glover 
et al. || find from their global analysis that it is impossible to fit both the CDF data for 
E T > 200 GeV and Deep Inelastic Scattering [DIS] data for x > 0.3. However, as noted 
in H the interpretation of large E T jet cross sections inherits uncertainties from the non- 
perturbative parton distribution and fragmentation functions. It has recently been reported 
Q that the apparent discrepancy between CDF data and theory may be explained by the 
uncertainties resulting from the non-perturbative parton distribution, in particular in the 
gluon distribution at large x. These authors have also performed NLO QCD global analysis 
including the CDF data and conclude that high E T can be explained in terms of a modified 
gluon PDF. However, we note that Lai et al. |3| use more parameters to describe the input 
gluon distribution than is usually done, whereas Glover et al. || assume that the gluon 
distribution has a canonical behavior [i.e. goes as (1 — x) n ] at large x. 
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It is also tempting to look for a possible explanation for the CDF data in terms of new 
physics. The CDF group have reported on a model of presence of quark substructure |J. 
They have compared their data to leading order QCD calculation including compositeness 
and have used MRSDO' parton distribution. They find a good agreement between data and 
the compositeness model, for E T > 200 GeV, for a substructure scale of Ac = 1.6 TeV. 

Yet another possibility is that jet measurements at hadron colliders may be sensitive 
to quantum corrections due to virtual SUSY particles |6|-f|. The purpose of this note is to 
concentrate on this scenario. The layout and of this paper is as follows. In next section 
we discuss our calculation of the SUSY virtual threshold effects. The final section contains 
discussions and conclusions of our numerical results. 



II. VIRTUAL SUSY THRESHOLD EFFECTS 



We consider the SUSY one loop corrections to the process d<r(pp — > 2 jets). As is 
well-known, the 2 jets production cross section in proton anti-proton collisions is found by 
weighing the expressions for differential cross section of the subprocess, d&y = da(ij — > 
2 final partons), by the parton distribution functions, and integrating over the parton vari- 
ables, i.e. 

da(pp -> 2 jets) = f dx x f dx 2 [f i / p (x 1 , Q 2 ) fj/ p (x 2 , Q 2 )]da i:j (a s , s, i, u). (1) 

Here da^ represents the subprocess cross section at cm. energy square of s = Xix 2 s, where 
y/s is the cm. energy of the pp system. It is well-known that SUSY particles [gluinos, 
squarks] decrease or slow down the rate of fall of a s (/j,) for large scale \i. "Large" means far 
above the threshold. At the one loop level the evolution equation for a s {p) can be written 

as 
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In the SM, b 3 is given by 



63 = 11 - g"/ , 



(3) 



whereas in MSSM model one has 
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where rif [nf\ is the number of quark [squark] flavors that are active. The contribution '—2' 
is the gluino contribution [it is assumed that the gluino is active in this case]. 

With this in mind, we can see that the SUSY corrections to Eq. ([!]) can be broken into 
three parts, namely (a): SUSY corrections to the PDF, (b): SUSY corrections to the running 
of a s , (c): parton- level SUSY loop corrections to d&ij, excluding the corrections already 
included in (b) [i.e. the running of a s ]. And finally we have to (d): combine those three 
part, as in Eq. ([!]), by convoluting with PDF to get hadron-level SUSY interference effects. 
In this note we consider all these four parts [i.e. (a),(b),(c),(d)]. Previous works on SUSY 
corrections |6|-[8| have considered only part of them: In Ref. || only the issue of the effect of 
high-mass thresholds due to gluinos, squarks and other new heavy quarks on the evolution 
of a s was considered [i.e. (b)]. The corrections to a s were found to be appreciable, this in 
turn means a significant increase in the transverse momentum dependence of jet production 
at the Tevatron. However, as noted in Ref. |7[], these authors || do not include the effect 
of qqg Yukawa interactions, and hence one cannot take their estimates for the superpartners 
of ordinary matter as final. In Ref. [0] the effect of Yukawa couplings was included, and 
found that the CDF data cannot be explained by a mass threshold effect in the MSSM, as 
the calculated result is not only small but of the wrong sign, considered at the parton level 
below the threshold energy scale [i.e. (b),(c)]. In a similar but more detailed analysis, the 
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authors of Ref. ||, working in the context of MSSM, consider at around the threshold energy 
scale the virtual one- loop corrections to the parton-level subprocesses qq — > qq, qq — > q'cf, 
qq' — > qq', qq — > gg and qg — > qg, which are expected to dominate the large E T cross sections 
at the Tevatron energy [i.e. (b),(c)]. 

The purpose of this note is to give our results of incorporating the one-loop radiative 
corrections into the running of a s , the dressing-up of the parton distribution functions, and 
finally convoluting the relevant subprocess cross sections with the SUSY dressed-up PDF's 
[i.e. (a),(b),(c),(d)]. In the hadron colliders, like the Tevatron, what is measured is pp cross 
section, and not the individual subprocesses cross sections. So in order to determine the 
effect of subprocesses on the E T cross section one must perform a convolution of the cross 
section of each subprocess with the corresponding PDF's [i.e. (d)]. We find that in the 
process of convolution with the PDF's, the "dips and peaks" in the various subprocesses |5) 
are much reduced. 

We note that, also as pointed out in ||, one should take into account sparticles effects 
on the parton structure functions at energies sufficiently far above the threshold, and 
can ignore this effect around the threshold region [i.e. (a)]. We test the validity of this 
statement and find it to be true from our numerical work, as will be shown in Fig. 2. We 
have considered the combined evolution equations for a s , q v 
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q NS (x,Q 2 ) 



a s (Q 2 ) 



P NS ®q NS (x,Q 2 ) . 



(6) 



dlnQ 2 



The second equation becomes, with the definition of p = xq 




(7) 
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The increase of a s (Q 2 ) results in the decrease of the PDF compared to the SM when we make 
the evolution. The qualitative reason is rather simple: stronger a s (Q 2 ) as Q 2 increases imply 
that the gluon radiation from the initial quarks are enhanced, and the PDF evolution yields 
the larger gluon or sea-quark densities at the smaller x region, and therefore the valence quark 
distribution q v (x, Q 2 ) decreases at large x, as a s (Q 2 ) increases, and vice versa. We recall that 
for large x the valence contribution dominates (e.g. q v ■ q v : q v ■ g : g ■ g ~ 0.65 : 0.3 : 0.05 at 
x ~ 0.3 at Tevatron energies), we ignored the SUSY evolution of sea-quark or gluon for large 
E T of CDF. It turns out, as shown in Fig. 2 and claimed in Ref. ||, the PDF with sparticle 
effects at around the threshold region deviates much less than 0.1% from its SM predictions, 
which justifies our assumption. And we can even totally ignore the SUSY corrections to 
PDF [i.e. (a)] all together for investigations below or around the threshold region. 



III. NUMERICAL RESULTS AND CONCLUSIONS 



For our numerical calculation, we implement various lower bounds on squarks and 



gluinos depending on parameters in the MSSM. For example, DO group |L0| searched for the 
events with large missing E T with three or more jets, observed no such events above the 
level expected in the SM. This puts some limits on the squark and gluino masses assuming 
the short-lived gluinos: 

m-g > 144 GeV for m~ q = oo , (8) 
or m g = nig > 212 GeV . (9) 

CDF group is currently analyzing their data, with their preliminary data being similar to 
the DO results with slight increase in sparticles' mass bound. In the subsequent numerical 



7 



analysis, we choose three sets of (nig, Trig), which we shall refer to as Case I, II and III, 
respectively, 

(m g , rriq) = (220 GeV, 220 GeV), (150 GeV, oo), (150 GeV, 150 GeV) . (10) 

The Case III is only of academic interest if the limit given in Eq. (|9]) is valid in reality. 

We turn now to our numerical result. The solid curve in Fig. 1 represents the one-loop 
evolution of a s (fi) versus the renormalization scale fi for the SM, assuming a top-quark mass 
of 175 GeV, and Aq C d = 0.2 GeV. The long-dashed curve represents the results for the 
MSSM with the five squarks and gluinos being degenerate with a common mass of 220 GeV 
[we shall refer to this as Case I]. As can be seen from Fig. 1, Case I starts to deviate from 
the SM results at fi = 440 GeV, as expected. The deviation in the value of a s away from the 
SM is maximum [for our range of at fi — 700 GeV and roughly on the order of 5%. This 
value is to be compared to the 12% reported in 0. The dotted curve in Fig. 1 describes 
our result for the MSSM with gluino mass of 150 GeV, and the squarks are assumed to be 
decoupled, we shall refer to this as Case II. Case II starts to be different from the SM at 
/i = 300 GeV. The maximum deviation for Case II occurs at around 700 GeV, where it also 
approaches the curve for Case I. 

The deviation from the SM of the one-loop evolution for the parton distribution func- 
tion [( sm" ) Uv (Q)] versus ^e factorization scale Q for Cases I and II are shown in Fig. 2. 
Case I deviates from the SM at Q = 440 GeV attaining a maximum change of —0.05% at 
700 GeV. The magnitude of deviation is much smaller compared to the results of Case I 
displayed in Fig. 1. The sign of corrections is also opposite to those for Case I in Fig. 1. 
This indicates as that a s increases [as for e.g. it does for MSSM] u v (Q) decreases. This is 
natural to expect since for larger a s there is more probability of gluon radiation, which in 
turn implies larger gluon and sea densities at small x, which eat away at the u v contribution 
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thus resulting in a reduction in u v . Turning to Case II, we see from Fig. 2 that it deviates 

from the SM at Q — 300 GeV, reaching a maximum value of only —0.1% at 700 GeV. 

Fig. 3 exhibits our results for [( MSS ^" SM )d~~] versus p T . for both Cases I and III [Case 

T j J 

III refers the MSSM with the five squarks and gluinos being degenerate with a common mass 
of 150 GeV, i.e. the last set in Eq. (|10|)1. A rapidity cut of 0.7 is used [\r)\ < 0.7]. A maximum 
difference of —1% from the SM is found for the both Cases. The question that naturally 
arises is, that why the "dips and peaks" [which are at the level of 5%-6% || ] at the parton 
level, are reduced. There are two factors contributing to the reduction in the "dips and 
peaks" in the process of going from the parton to the process level. Part of this reduction 
comes from the process of convolution of the various subprocesses with PDF's [as already 
remarked to in the previous section]. The other piece of this reduction in "dips and peaks" 
by going from subprocess to process level is due to t-channel "dilution" effect. This can be 
shown by simply not including the t-channel subprocesses' contribution. When this is done 
[see Fig. 4a and 4b] reduction in the "peaks and dips "is not so large. 

In Fig. 4(a) we exhibit our results for [( MSS ^~ SM )^-] versus p T . for both Cases I 
and III, but this time only including the subprocesses qq — > q'q' and qq — > gg. These 
two are the subprocesses with the prominent threshold structures from vacuum polarization 
interferences through s-channel exchange diagrams. A rapidity cut of 0.7 is used [\rj\ < 0.7]. 
The deviation from the SM for Case I varies between 2% to —4%, whereas as in Case II the 
variation is similar size to Case I at different energy scales related to the SUSY threshold 
effects. In Fig. 4(b) we exhibit our results for [( MSS S ^~ SM )^] versus M js for both the Cases 
I and III, including the subprocesses qq — > q'q' and qq — > gg. The percentage variation is 
almost the same as in Fig. 4(a). 

Near the virtual (or direct) SUSY threshold, the Coulomb interactions between the 
gluinos [or squarks] with a very small velocity v = Jl — Arn-fQ 2 make the QCD effects 



9 



depending on the parameter a s /v, rather than on a s . The result of summing up the Coulomb 



exchanges reduces to multiplying the Coulomb factor a s /v |jTT|, which gives a wide resonance 
structure around Q « 2rrig for the short-lived gluino pairs. If the life time of gluino is long 
enough, then narrow bound states of gluino pair may appear near the threshold region. Since 
the properties on gluino pair bound states are not known and largely model dependent, we 
can guess only qualitative nature of the threshold region, Q w 2m g (or Q « 2m^) , as shown 
in Fig. 3 and 4. 

In summary, in the MSSM the E T distributions does not differ very much from those of 
the SM except for the possible threshold effects (~ 1%) through loop corrections. In actual 
experiments, the jet resolution will in general smear out any narrow resonance structures 
(which may be the case for the long-lived gluinos), leading to broad resonance structure, and 
therefore it looks impossible to detect SUSY particles through this kind of indirect virtual 
threshold effect. As is previously explained, it has been reported Q that the apparent 
discrepancy between CDF data and theory may be explained by the uncertainties resulting 
from the non-perturbative parton distribution, in particular in the gluon distribution. Our 
main conclusion seems to be on the right track in view of this global analysis of parton 
distribution function. 
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FIGURES 

Fig. 1 The one- loop evolution of a s (fi) versus /i is given for SM [solid line] and SUSY cases. 

Fig. 2 Percentage deviation from the SM due to SUSY contribution of PDF of the valence 
quark versus factorization scale Q for typical momentum fraction x = 0.3. 

Fig. 3 Deviation from the SM due to SUSY contribution to 4 s - versus p T . 

Fig. 4a Deviation from the SM due to SUSY contribution to 4 s — versus r>„ , when only the 

ap T , 1 j 
3 

subprocesses qq — > q'q', and qq — > gg are included. 
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Fig. 4b Deviation from the SM due to SUSY contribution to versus Mjj, when only the 
subprocesses qq — > q'q', and qq — > gg are included. 
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